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Lattice QCD is used to calculate the potential between a static quark and anti-quark in the presence 
of a finite density of 7r + 's. Correlation functions of multiple 7r + 's are used in conjunction with 
Wilson-loop correlators to determine the difference between the QQ-potential in free-space and in 
the presence of a pion condensate. The modifications to the potential are found to have significant 
dependence on the QQ separation over the range r < 1 fm explored in this work. Our results are 
consistent with the pion-condensate behaving as a (nonlinear) chromo-dielectric. 



Introduction: Flavor dynamics plays a central role in the 
study of heavy-ion collisions and, ultimately, in unravel- 
ing the phase structure of matter at finite temperature 
and density. The production, evolution and detection of 
the heavy flavors, charm and bottom, in such collisions 
has long been known to be sensitive to the presence of 
deconfined phases of Quantum Chromo dynamics (QCD), 
such as the naive quarkjduon plasma or "the perfect 
fluid" seen at RHIC [H, [J. Color charges are screened 
in a deconfined phase, and the production of a given 
state of quarkonium is expected to be suppressed when 
the screening length becomes comparable to, or smaller 
than, its size 3] (for a recent discussion see Ref. 0]). 
However, quantitative calculations of " .//^-suppression" 
and more generally "quarkonium-suppression" in heavy- 
ion collisions are difficult due to the enormous complexity 
of such collisions. 

Understanding the behavior of a static quark, Q, and 
static anti-quark, Q, pair produced in a heavy-ion colli- 
sion requires knowledge of its dynamics in both the de- 
confined and hadronic phases. In the hadronic phase, the 
potential between a static quark and static anti-quark 
(denoted as the QQ potential) can be screened by the 
hadrons that participate in the collision. Differences be- 
tween this "hadronic medium effect" and the effect of a 
deconfined phase are tell-tale signatures of a deconfining 
phase transition in heavy-ion collisions. As a step to im- 
proving the description of QQ transport in the hadronic 
phase we present the results of a Lattice QCD calcula- 
tion of the modifications of the QQ potential and force 
resulting from the presence of a condensate of charged 
pions (either all 7r + 's or all 7r~'s) in the absence of elec- 
tromagnetism. The QQ-potential is found to be lowered 
in the condensate, with a non-trivial dependence upon 
the QQ separation, leading to a reduced QQ force. 

It is useful to consider the behavior of the in-medium 
QQ potential and force in the limiting cases where the 



separation, R, is much smaller than, or much greater 
than, the scale of chiral symmetry breaking, A x , that 
typifies hadronic interactions. In the R A x -C 1 limit, 
it is appropriate to construct an effective field theory 
(EFT) describing the interactions between the QQ with 
the gluon fields. The coefficients in this EFT are deter- 
mined via a multi-pole expansion of the matrix elements 
calculated in QCD. The interactions between the QQ- 
state and the hadronic background are thereby factorized 
into short-distance parts encapsulated in the coefficients 
of the EFT, and low-energy matrix elements of local op- 
erators composed of quark and gluon fields (explicitly 
independent of R). This method [5| has been used to 
calculate the binding energy of quarkonium to infinite 
nuclear matter [6[, and quarkonium scattering from nu- 
cleons Q. The interactions of a spatially-averaged and 
orientation-averaged Wilson-loop with quark and gluon 
fields, at leading order in the strong-coupling and deriva- 
tive expansion, are described by an effective Lagrange 
density of the form 



C = R 3 S^S G^Gl [ asT + c 2 v»v u } 



(1) 



where S is the operator that annihilates the QQ-pair, 
G a (3 is the gluon field strength tensor, v a is the QQ- 
pair four-velocity. The renormalization-scale dependent 
coefficients Cj that appear in eq. |T]) are dimensionless and 
are expected to be of order unity by naive dimensional 
analysis. The Lagrange density in eq. Jl} makes explicit 
the separation dependence of the in-medium component 
of the potential (oc R 3 ) and force (ex R 2 ) for R A x <C 1. 

At extremely large distances, the ground state of the 
system is a heavy meson-anti-meson pair in the pion con- 
densate interacting through the exchange of hadrons de- 
scribed by an effective field theory analogous to that used 
to describe the interactions of nucleons. In QCD, the 
force between these heavy mesons is Yukawa-like with a 
mass-scale set by the pion mass. 
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At intermediate distances and in vacuum, the ground- 
state of the system has the QQ potential increasing with 
separation, and the gluonic field configuration between 
the QQ is tending toward a flux-tube with constant force 
between the QQ pair (for recent lattice calculations see 
Ref. [8j). It is not obvious how the presence of a pion- 
condensate will modify the interactions between the QQ 
pair. If instead, the system under consideration was a 
collection of neutral hadrons in the presence of an electric 
field, the energy-shift of the system would depend upon 
the electric-polarizabilities of the hadron and the volume 
of the electric field. An overall reduction in the QQ force 
would result, encapsulated in the dielectric function of 
the medium. 

The Lattice QCD Calculation: The ground-state energy 
of a system comprised of a QQ-pair separated by a dis- 
tance R and n-ir + 's can be extracted from the correlation 
functions 

c„(W) = (o|[XIx»+(3c,*)xt + (o,* ff )] n |o) , 

X 

C w {R,t w ,t) = (o| W(y + r,t;y,^)|o) , 

c n , w (R,U,t w ,t) = ^o| [^X 7r +(x,t)xj r +(0,< lr 





AWp-0 m AWp-1 "* AWp-2 * Nuyp-4 



x W(y + r,i;y,0 

y,\r\=R 



(2) 



where Xk+( x ) — ua {x)lbd a {x) is a pseudoscalar inter- 
polating operator for the 7r + (a is a color index), and 
W (y, to; y + r, t) is the Wilson-loop operator formed 
from products of gauge-links joining the vertices at 
(y,*o), (y+J",*o), (y + r,t) and (y,t). At large times, 
the static QQ potential in vacuum is extracted using 



t-u, 



Cw(R> t Wl t w + 1) 



(3) 



It is useful to define the ratio of the three correlation 
functions in eq. @, 



Gn,W \R) ^7r , t w , t) 



Cn(tn,t) Cw{R,t w ,t) 



(4) 



from which the in-medium modification to the potential 
can be extracted, for t w 3> t n (ensuring that the pion 
system is in its ground-state), 



(log 



G 



5V(R,n) ,(5) 



where the (•••)*„ denotes an average over a number of 
initial time-slices for the Wilson-loop, i™ m to i™ ax . 

We have computed the correlators in eq. ((2]) in mixed- 
action lattice QCD, using domain- wall valence quark 
propagators from a Gaussian smeared-source on rooted- 
staggered MILC gauge-configurations (see Refs. 0, [l(| 
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FIG. 1: The QQ force in vacuum determined with different 
levels of HYP-smearing. Inset is the QQ potential with these 
same smearings, each normalized to V(Ro) = at Rq — 5b. 



for details). Here we focus on calculations on an ensem- 
ble of 1001 lattices with a pion mass of ~ 320 MeV, 
and a lattice spacing of b = 0.087(1) fm with dimen- 
sion 28 3 x 96 giving a spatial dimension of ~ 2.5 fm. 
Propagators were calculated with both periodic and anti- 
periodic boundary conditions in the time-direction and 
combined to effectively double the length of the time- 
direction leading to long plateaus in effective energy plots 
of the mesonic correlators. The light quark propagators 
were computed after the gauge-field had undergone a sin- 
gle level of HYP-smearing [11[ , while Wilson-loops were 
calculated with Ahyp = 0,1,2 and 4 levels of HYP- 
smearing of the gauge-field. Further, the spatial links 
of the Wilson-loops were APE-smeared [12, EH in the 
transverse-spatial directions to optimize the signals for 
Cw{R, t w , t). The increasing levels of HYP-smearing re- 
sult in improved signal-to-noise ratios and enable the po- 
tential, and in-medium modifications to the potential, to 
be determined over a range of QQ separations. The cor- 
relation functions C„ (t^ , t) have been previously calcu- 
lated on these lattices, enabling a study of the properties 
of the pion and kaon condensates, and the three-meson 
interactions 3, HI, 16j |. 



Correlated fits are performed to the effective energies 
(or energy-differences), e.g. eq. ([5]), for each QQ sepa- 
ration and number of pions in the volume. Separately, 
jackknife and bootstrap procedures were used to gener- 
ate the covariance matrix over the given fitting interval 
of time-slices, and correlated x 2 - m i nm h za tion was per- 
formed to extract the energy and its associated statisti- 
cal uncertainty. A systematic uncertainty is determined 
by a comparison of the various fit procedures and various 
fitting ranges (including the choice of i™ m ^ max ). 

The QQ-potential in vacuum is determined from Cw, 
and the QQ-force is determined by finite-differences of 
the potential. Our results for these two quantities are 
shown in fig. [T] By comparing the force in vacuum calcu- 
lated with different levels of HYP-smearing, we ascertain 
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FIG. 2: Representative effective energy plots for the in- 
medium contributions to the potential. The horizontal band 
corresponds to the combined statistical and systematic uncer- 
tainty of the respective fit. 
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FIG. 4: The in-medium contribution to the QQ potential 
at R ~ 0.24 fm as a function of the number of pions. The 
uncertainties are as described in fig. [3] The isospin chemical 
potential of the system is shown on the upper axis flfil ]. 



the separations at which the force (and potential) cease 
to be significantly contaminated by the smearing proce- 
dure. We conclude from fig. [1] that the potential and 
force calculated at separations R > bN^yp are close to 
the result of an un-smeared calculation. In our analysis, 
we only use calculations that satisfy this criterion. Three 
representative effective energy plots associated with the 
in-medium contributions to the potential are shown in 
fig. O Typically, signals become noisier as either R or n 
increases, restricting the present analysis to R < 1 fm. 

The in-medium contributions to the QQ-potential, 
8V{R 1 n), resulting from both a single pion in the lat- 
tice volume (a number density of po ~ 1/(2.5 fm) 3 
0.0 f i 1 fm -3 ), and from five pions in the lattice volume 
(p = 5po = 0.32 fm -3 ), are shown in fig. [3] For p < 7po 
they are found to be linear in the density within the un- 
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FIG. 3: The in-medium contributions to the QQ potential 
at a pion densities of po and 5po- The inner uncertainty as- 
sociated with each point is statistical, while the outer is the 
statistical and systematic uncertainties combined in quadra- 
ture. Different HYP-smearings are offset for clarity. 



certainties of the calculations, as can be seen for a repre- 
sentative R in fig. [4] We note that in-medium contribu- 
tions to the potential as small as 5V(b, 1) ~ 100 keV 
have been determined. For the system containing a 
single pion in the lattice volume, the energy-shift can 
be directly related to the scattering phase-shift using 
Luscher's method |17j. Therefore, the 8V(R,1) shown 



in fig. [3] can be used to determine the scattering length 
associated with a QQ-pair of fixed separation and a pion. 

The effects of the medium on the radial QQ force, 
SF(R,n), are determined from the effective energy- 
differences derived from finite differences of 5V(R, n) 
with respect to R. The modifications to the force at 
densities po and 5po are shown in fig. [5l The QQ force is 
seen to be reduced by an approximately i?-independent 
amount over the separations and pion densities we have 
been able to explore. 

The in-medium potential and force have also been cal- 



12 - 
10 1 

8 1 
6 1 
4 1 
2 1 
1 



n=5 




ifi yv HYP =o i§i yv HYP =l ifi A^fYp=2 i^i N H 



o.l 



0.2 



0.3 0.4 
R [fm] 



0.5 



FIG. 5: The in-medium contribution to the radial QQ force, 
5F(R,n), at a pion number density of po and 5po. The un- 
certainties are as described in the caption of fig. [3] 



4 



culated on MILC lattices of the same spatial volume 
as the configurations used here, at four different quark 
masses, but at a coarser lattice spacing (b ~ 0.125 fm). 
We find the results at the corresponding pion mass are 
consistent, suggesting that lattice discretization errors 
are not large. However, the uncertainties in the medium 
modifications to the QQ potential and force are some- 
what larger than those on the current ensemble. A 
mild variation with the pion mass was observed, how- 
ever, more precise calculations are required in order to 
quantify this dependence. 

Discussion: In-medium effects play an important role in 
the diagnostics used to explore new phases of matter in 
heavy ion collisions, and more generally emerge as useful 
quantities in the context of mean-field constructions in 
many-body systems. We have performed the first QCD 
calculation (without electromagnetism) of modifications 
to the QQ potential and force in a hadronic medium by 
calculating the energy of a Wilson-loop in the presence 
of a condensate of charged pions. The attractive QQ in- 
teraction is found to be reduced by the hadronic medium 
over the range of separations we were able to explore, 
R < 1 fm. This is a first step toward a more systematic 
exploration of in-medium effects with Lattice QCD, with 
the ultimate goal of looking for in-medium modifications 
of hadronic observables in backgrounds of baryons. Such 
calculations will require significantly more computational 
resources than are currently available, and precise calcu- 
lations of multi-baryon systems are a prerequisite. 

At p , our calculations are analogous to those of J/ip-n 
scattering lengths (with the J/ip replaced by a Wilson- 
loop) which have been performed previously in QCD [181 ] 
and quenched QCD [191] , At higher densities, the calcula- 
tions involve multi-pion backgrounds and are the first of 
their kind. A non-zero three-pion interaction is required 
to describe the volume-dependence of the energy-levels of 
n > 2 7r + 's in these volumes 3, HI 3], however, within 
the uncertainties of our calculations, the multi-pion in- 
teractions with the QQ-pair are found to be consistent 
with zero over the range of QQ separations we have ex- 
plored. It is important to refine this work by performing 
higher-statistics calculations in order to determine the 
multi-pion interactions with the QQ-pair. In addition, 
with a corresponding calculation of the pionic matrix el- 
ements of the gluonic operators in eq. (fTJ), the coefficients, 
Cj, in eq. |T]) could be determined. 

From our analysis we observe that the medium- 
modification of the force is independent of the QQ sepa- 
ration in the region where the force in vacuum is becom- 
ing constant. This is consistent with the pion-condensate 
behaving as a dielectric in the volume of the color flux- 
tube between the QQ-pair. It implies that the pion and, 
collectively, the condensate has a chromo-susceptibility, 
which is expected to be highly non-linear in the gluon- 
field strength. At small separations, the modification to 
the force appears consistent with the behavior expected 



from eq. ((T|), but calculations at smaller lattice spacing 
are required to confirm this. 
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